Effects of the Tensor Force on the Multipole Response in Finite Nuclei 
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We present a thorough analysis of the effects of the tensor interaction on the multipole response 
of magic nuclei, using the fully self-consistent Random Phase Approximation (RPA) model with 
Skyrme interactions. We disentangle the modifications to the static mean field induced by the 
tensor terms, and the specific features of the residual particle-hole (p-h) tensor interaction, for 
quadrupole (2"'"), octupole (3~), and also magnetic dipole (1^) responses. It is pointed out that the 
tensor force has a larger effect on the magnetic dipole states than on the natural parity states 2+ 
and 3~ , especially at the mean field level. Perspectives for a better assessment of the tensor force 
parameters are eventually discussed. 

D ' PACS numbers: 21.60.Jz, 21.65.Ef, 24.30.Cz, 24.30.Gd 

(S| ; I- INTRODUCTION 

The nuclear effective interactions like the zero-range Skyrme forces have been quite successful in many respects 
since they have been introduced, mainly because of continuous efforts to apply them to the study of different nuclear 
phenomena, and to improve their form and their parameters. The Skyrme forces contain typically 10 free parameters 
^ , which are fitted to reproduce empirical bulk properties of uniform nuclear matter, and masses and radii of a few 
(— ( ' magic nuclei. These forces, or at least several of the many existing parameter sets, describe in a reasonable way the 

' global features of the ground-state properties along the nuclear chart (i.e., binding energies, radii, deformations). In 

some cases one needs to input more information in the parameter fitting: the Lyon forces [l| have been for instance 
determined by requiring a reproduction of the equation of state of neutron matter emerging from realistic calculations. 
1 Properties of excited states (vibrations, rotations) have been studied using the Skyrme parameter sets, allowing a 
■ large amount of physical insight. For a review on mean-field calculations the reader can consult Ref. Few ideas 
t;^ , about selecting the Skyrme forces which have reasonable overall performances can be found in Last but not least, 
in the quest for a universal local Energy Density Functional (EDF) for nuclei, the Skyrme framework is often used as 
, a starting point. 

' Of course, there are drawbacks and limitations of the mean-field approaches based on effective forces. This issue 
^\ \ should be discussed together with the answer to the question whether we need to generalize the existing Skyrme 
• parametrizations. In this spirit, since most of the Skyrme sets which have been widely used are purely central, many 
groups have recently devoted much attention to the role played by the zero-range tensor terms that can be added to 
the standard Skyrme terms (see Refs. [1, 0, H, [1, E^j EH)- This blooming of theoretical studies has followed 
' the claim by the authors of Ref. that the tensor force is crucial for the understanding of the evolution of the 
single-particle energies in exotic nuclei. 

Certainly, the mentioned studies have elucidated that the tensor force does play a role. At the same time, the 
single-particle states are not the right observable on which to constrain an effective force or an effective functional, 
since for these states other correlations beyond mean-field are expected to manifest themselves. 

It is certainly quite timely to analyze in detail how the tensor terms added to a standard Skyrme force, or the new 
parameter sets introduced in |8i| (fitted with the tensor included), behave as far as the excited states are concerned. 
The study of nuclear collective vibrations, within self-consistent Random Phase Approximation (RPA) on top of 
Hartree-Fock (HF), is a well-defined framework which has allowed testing the effective Skyrme sets for many years. 
Recently, we have developed fully self-consistent RPA JJ5.], and quasi-particle RPA (QRPA) 14], schemes. The full 
self-consistency, in the nonrelativistic framework, is essential to respect theorems concerning sum rules associated 
with appropriate operators. In this paper, we extend these schemes by including the tensor terms both in HF and in 
RPA as terms of the residual interaction, and we analyze the results in the case of the multipole response of magic 
nuclei. 

Our work is the first one which attempts a detailed explanation of the effects induced by the tensor force on the 
response of finite nuclei within the Skyrme framework. However, the fact that the subject is of interest at present 
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is testified by recent papers that are devoted to related topics. V. De Donno et al. have studied the low- 



lying magnetic excitations in several nuclei, either within a phenomenological approach or by using the Gogny force. 
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The authors of [16[ have investigated self-consistently, using the Skyrme forces and the zero-range tensor terms, the 
response of uniform matter. 

The outhne of our paper is the following. Since Skyrme- RPA is a well-known framework we will survey very briefly 
our formalism in Sec. [TTl More emphasis is given to the discussion of the results, for quadrupole, octupole and 
magnetic dipole states, in Sec. IIIII An overall conclusion, together with some perspectives for future work on the 
tensor force and on possible constraints on its parameters, is given in Sec. IIVI The main technical part of this 
work, namely the evaluation of the p-h matrix elements of the tensor interaction, is described in Appendix |X1 while 
Appendix |B] describes and approximate yet quite useful separable approximation for the tensor force that is used as 
a guideline to understand some features of our numerical results. 



II. FORMALISM 



As mentioned in the Introduction, Skyrme-RPA theory is well known since many years, especially in its matrix 
formulation. We have been using for several applications a scheme which is fully self-consistent, that is, where there 
is no approximation in the residual interaction since all its terms are taken into account including the two-body spin- 
orbit and Coulomb. In our scheme we use box boundary conditions. In other words, the continuum is discretized. In 
the present case the box radius is set at 15 fm for the Ca isotopes and at 20 fm for ^"^Pb. 

After solving the HF equations in coordinate space and determining, together with the mean field, the unoccupied 
states as eigenstates of this mean field in the mentioned large box, we build up a model space of p-h configurations 
with given multipolarity J'^ and we write and solve the RPA matrix equation in that space. The model space includes 
the configurations built up with all hole states, and with the particle states labelled by rip, Ip and jp (number of radial 
nodes, spatial and total angular momentum); Ip and jp take all possible values which are allowed by selection rules, 
while Hp varies between nmax + 1 (T^-max is the number of nodes of the last occupied state) and nmax + 1 + An. 

The value of An is chosen large enough, so that the isoscalar (IS) EWSR for the multipoles 2+, 3~ exhaust 
practically all (~ 97-99%) the double commutator (DC) value. In particular, we use An=10 unoccupied shells in the 
case of the 2+ and 1+ calculations in ^"Ca and ''^Ca, An=9 unoccupied shells in the case of the 3~ calculations in 
the same nuclei, and An=8 unoccupied shells for ^°®Pb. 

The eigenvalues £"„ and eigenvectors |n) of the RPA equations allow calculating the response function to different 
operators. In particular we are going in what follows to show figures with the strength function associated with the 
operator F, namely 

^(i?) = ^|(n|F|0)|2<5(£;-i?„), (1) 

n 

where the sum spans the whole RPA spectrum. In the case of discretized RPA one has in reality sharp peaks but 
in the figures we shall display S{E) smoothed out using Lorentzian functions (having 1 MeV width). We shall use 
isoscalar quadrupole and octupole operators, that is, 

A 

Fl = Y. V2rTT r^YLMih) (2) 

i=l 

with L=2 and 3, respectively. We focus our study on the 1+ spin-flip states as well. In this case the following operators 
are used, namely 

A IS k 

^•^(^^) = E ^VJ(2^+1) <yoo{r.) (3) 
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in the isoscalar case, and 



Fj{IV) = <Yooin)r:- (4) 



in the isovector case. We shall consider the case J—1. In the above equations the nuclear magneton ii^ = 
appears together with the quantities g^'^ ^ jigs + 9s - 1)=0.19, and g^^ ^ jigs - 9s + 1)=-2.10. 

We note that the tensor contribution must not alter the DC value for the EWSRs in the case of the operators 
for all values of L. It does change, instead, the value of the EWSRs in the case of the operators ^ and (jH) (cf., e.g., 
Ref. flT^). 
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Two sets of Skyrme parameters are employed in the calculations presented below. On the one hand, we have chosen 
the set SLy5 Jj and supplemented it with tensor terms characterized by the values chosen in and used in as 
well. The analysis of these results will elucidate qualitatively the role played by the tensor terms, but in order to 
compare with the case of a force having all parameters fitted on equal footing (including those of the tensor terms) 
we have performed also calculation with the set T44 introduced in Ref. t^]. For the reader's convenience, all the 
parameters are shown in Table [H 

III. RESULTS 
A. Quadrupole response 

Fig. [T] displays in the upper panel the IS response of the nucleus ''"Ca. For a I ■ s-closed nucleus (that is, a nucleus 
where both spin-orbit partner states are occupied), and a natural-parity excitation spectrum, we do not expect that 
the tensor interaction plays any significant role. In particular, for the ground state, the changes in the spin-orbit 
splittings due to the tensor terms have been written in previous works and they depend on the so-called spin-orbit 
densities defined as 

Mr) = ^ E(2J'» + + 1) - ^'(^^ + 1) - il^'W' (5) 

i 

where the isospin quantum number q = 0(1) labels neutrons (protons) and the index i — n,l,j runs over all states 
having the given q. Ri{r) = is the radial part of the wavefunction. So, our expectation is obviously based, in 
the unperturbed case, on the fact that Jq is negligibly small (both for neutrons and for protons) in the case at hand. 
In fact, the microscopic calculations show that the results are exactly the same, with and without the tensor force. 

For comparison, we have computed the same IS response in the nucleus "^^Ca and the results are shown in the 
middle panel of Fig. [1] The response in the giant resonance region is not affected strongly by the inclusion of the 
tensor force. Although some changes are found in the high-energy part of the strength function, the main impact 
of the tensor terms is visible in the states below the isoscalar giant quadrupole resonance (ISGQR). For the sake of 
clarity, the properties of the lowest 2+ states are shown in Table [Til We have analyzed the changes of these properties 
induced by the tensor force. With the interaction SLy5, and the tensor parameters T and U already used in Refs. 
[1, the spin-orbit splittings are increased. Normally, at least for well-bound states, the splittings should increase 
more for larger values of I. Consequently, the low- lying 2+ state in ^*Ca, which is mainly due to the neutron f7/2 — > 
P3/2 transition, is pushed upward in the unperturbed response by the tensor force. We denote this shift by AEnp. 
The effect of the residual tensor force Vtcnsor included in RPA can be estimated by means of 

AEjiPA « AEuF + (Vtcnsor) , (6) 

where AE'rpa indicates the difference between the RPA result with and without the tensor force, and () means that 
we extract here an average value of the residual force. In our SLy5 calculation, the value of A_Ehf is 1-41 MeV and 
from the shift of the RPA peak (0.83 MeV) we extract (Vtcnsor)=-0.58 MeV. 

The fact that the residual force is attractive, albeit not large, can be understood using the argument developed in 
Appendix [B1 Using a separable approximation for the tensor interaction, it is shown that both the tensor-even and 
tensor-odd terms produce indeed attraction in the present case, where T is positive and U is negative [cf. Eqs. (IB9P 
and ((BT3l) ]. 

Our understanding needs to be confirmed by another example. In the lower panel of Fig. [Jl the IS response for 
the nucleus ^"^Pb is displayed. There is essentially no effect of the tensor force in the giant resonance region. The 
low-lying 2+ state is moved upward in the calculation with tensor by A£'rpa=0.17 MeV. However, if we look at the 
centroid of the low-energy unperturbed response, the shift due to the tensor force is Ai5HF=0.60 MeV. Then, from 
Eq. ([6]) we extract a value of (Vtonsor)=-0.43 MeV, that is, consistent with what extracted in the case of "^^Ca and 
with the qualitative argument of Appendix [Bl 

The inclusion of the tensor force on top of SLy5, with the parameters used in our previous works (gI. [lOj. does not 
improve the agreement with experimental values. It is desirable to check the behavior of a different interaction, and 
we have chosen one of those which have been fitted including the tensor terms in Ref. [8], that is, the set named T44. 

In this case, not so much emphasis should be put on the discussion of the effect of the tensor terms separately from 
that of other Skyrme parameters, since all parameters have been fitted together. We just briefly point out that the 
tensor residual interaction is almost negligible in this case. In fact, for the low- lying 2+ state of ^^Ca calculated with 
T44, the values of the quantities defined in Eq. (HD are Ai;HF=-0.13 MeV, Ai;RPA=-0.18 MeV and (Vte„sor)=-0.05 
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MeV. The same values are A£;hf=-0.54 MeV, A£;rpa=-0.55 MeV and (Vtonsor)=-0.01 MeV in the case of ^ospb. We 
should also mention that in the case of infinite matter [iQ] the effect of the tensor terms has been found to be small, 
in the case of the Tij parameter sets of Ref. Q, looking at the S — response both with and without the tensor. 

In conclusion, the behavior of the tensor interaction in the quadrupole response can be fairly well understood. Its 
effects are visible only for the low-lying states (not for giant resonances) and are mainly related to the mean field 
since the residual tensor interaction matrix elements are either tiny or anyway not very large. We end this Section 
by adding that the isovector response has also been computed, although it is not displayed in our figures. We have in 
fact observed that the centroid energy of the (high-lying) isovector strength may change by few hundreds of keV and 
only in few cases by about w 1 MeV when the tensor force is included. 

B. Octupole response 

In a similar style as for the quadrupole case, the IS octupole responses of the nuclei ""^Ca, '*^Ca and ^"^Pb are 
displayed in the upper, middle and lower panels of Fig. [21 respectively. Also in the present case, we do not discuss 
the IV strength (for which a comment similar to the one made at the end of the previous subsection applies) but only 
the IS strength. 

The detailed behavior of the low- lying states is shown in Table IIIII In the case of ^"^Ca, as we have already 
commented in the previous subsection, there is no effect of the tensor terms at the mean field level since Jp and Jn 
are negligibly small. In the case of SLy5 plus the T and U values of [1, [13], the value of (Vtcnsor) can be therefore 
deduced from the shift of the RPA energy of the lowest 3^ state, and it is -0.76 MeV, not much different from the 
values obtained in the case of the low- lying 2+ (cf. the previous subsection). 

In the case of **Ca, there is no systematic behavior in the shifts of the unperturbed p-h transitions when the tensor 
terms are included. In this nucleus Jp is negligibly small while Jn is positive due to the neutron contribution. 
With the force SLy5, and the parameters T and U of [1, the proton spin-orbit splittings are reduced, whereas the 
neutron ones are enlarged. The 3j~ state, which appears without tensor at 4.78 MeV as it is shown in Table [1111 has 
the proton Si/2 h/2 E^iid d3/2 — > h /2 main components, and both are pushed upward in energy when the tensor 
is included. This explains the large positive shift of this RPA state, despite the (not strong) attractive contribution 
provided by the residual interaction. The same p-h configurations are also the main ones which contribute to the 
second RPA state at 5.80 MeV which is visible in Fig. [2l(in the curve corresponding to the calculation without tensor). 
However, in the case of other RPA states the situation is different. The third peak visible in Fig. [21 (once more, in 
the curve corresponding to the RPA calculation without tensor) lies at 9.23 MeV, and it receives contribution both 
from configurations which are pushed upward by the tensor inclusion, like the proton — > P3/2 and the neutron 
i-j/2 g9/2i a-iid from others which are pushed downward like the neutron d3/2 P3/2- So, the inclusion of the tensor 
cannot produce an effect which is understandable in simple terms. On top of this, we have verified that the RPA 
states calculated with the tensor included, contain somewhat different admixtures of p-h configurations with respect 
to those calculated without the tensor force. 

Also in ^"^Pb several configurations contribute to the collective low-lying 3^ state. This collectivity is enough, so 
that the effects coming from unperturbed configuration which are pushed upward or downward by the tensor terms, 
together with the effect of the tensor residual interaction, essentially cancel one another. Both from Fig. [2land Table 
IIIII one can see that the energies and B(E3) values in ^ospj^ affected significantly when tensor is included. 

The conclusion which can be reached, looking at the results obtained by using the SLy5 interaction and the tensor 
parameters from [l^l , is that the effects of the tensor force on the octupole states are rather dependent on the 
specific nucleus, that is, on the relevant shell model states. In any case, we do not observe huge effects. 

If we use the interaction T44, these qualitative conclusions remain valid, although the behavior of the tensor force 
is different in detail because of the contributions with different signs of T and C/, and their smaller absolute values. 
Actually, as in the quadrupole case, these differences do not manifest dramatically in the final results of the low-lying 
states in "'^Ca and ^"^Pb. In the case of ''"Ca, the gap in this nucleus turns out to be very small (about 3.3 MeV) as 
compared to the case of standard Skyrme forces which are fitted without tensor terms. In fact, this gap is, for example, 
5.51 MeV with SLy5. This explains the small energy of the low- lying 3~ in ^"Ca. In some cases this behavior of T44 
might lead to RPA instabilities, which, however, we have not observed by studying these low multipole responses in 
finite nuclei. 



C. Magnetic dipole response 

On quite general grounds, it can be expected that the effects of the tensor force are larger for spin (or for spin- 
isospin) states, both at the mean field level since unperturbed configurations are sensitive to the spin-orbit splittings. 
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and also as far as RPA correlations are concerned. Some of the authors of the present work have already analyzed 
the role of the tensor correlations in the case of charge-exchange, spin-isospin states like the Gamow- Teller (GT) and 
spin-quadrupole resonances [13] ■ Tensor correlations are quite strong in that case, and lead to a lowering of the main 
GT peak by about 2 MeV in ^"Zr and ^°^Pb, which is accompanied by a strong repulsive shift of a sizeable fraction 
of strength (pushed at energies higher than the giant resonance region, at 30 MeV and above). These effects are 
certainly more prominent than what has been discussed in the previous subsection for the non spin- flip transitions. 
Many earlier works in the literature have emphasized the role of the tensor force in the spin and spin-isospin channel, 
but not in a self-consistent framework like in the present case or in the case ofRef. H. The self-consistent calculation 
of the uniform matter response performed using the Skyrme force T44 in jl^ , confirms that larger and non-trivial 
effects from the tensor terms may be expected in the spin channel. With this background, we have made calculations 
for the magnetic dipole response in ^^Ca and ^°®Pb. 

The results we have obtained in the case of ^*Ca are displayed in Fig. O In the left panels, the strength functions 
associated with the SLy5 Skyrme set plus the T and U parameters of [6| are shown. The HF peak is associated with 
the neutron i-j/2 ^ h/2 configuration, which lies at 7.06 MeV without the tensor contribution, and is pushed at 10.68 
MeV by the inclusion of tensor terms. In the case of RPA, we deal with a somewhat undesired feature of the SLy5 
set, which is characterized by a positive value of the Landau parameter Go=1.14 and a negative value of Gq=-0.15. 
Repulsion in the spin-isoscalar (IS) channel and attraction in the spin-isovector (IV) channel is in contrast with the 
empirically accepted values of the Landau parameters, and, as already noticed in plj . constitutes an anomaly of 
the SLy5 set. In the case at hand, the residual interaction is essentially the J'^=\'^ diagonal matrix element of the 
mentioned neutron configuration, that is, it is the sum of the IS and IV part. Without tensor, from the fact that the 
RPA peak is at 9.28 MeV one can extract that the interaction is dominated by the IS contribution and is repulsive 
(its value being 2.22 MeV). Including the tensor interaction, the RPA peak moves to 12.31 MeV, so the shift due to 
tensor correlations is 3.03 MeV: the tensor residual interaction can then be extracted from Eq. ^ and it turns out 
to be -0.59 MeV. The fact that it is attractive is understandable in keeping with Eq. (jB18|) and from the negative 
sign of U . 

In the case of the T44 interaction, the values of the Landau parameters Go and Gg are, respectively, 0.40 and 0.06. 
The results associated with this parameter set are visible in the right panels of Fig. [31 and one can see that the 
effect of the tensor is quite small as far as both unperturbed mean-field and RPA correlations are concerned. This 
is related to the smallness of the parameter U [cf. Table [H Eq. ^ and the formulas of Appendix B]. In the case of 
the magnetic dipole state we have, therefore, found that different ways of including the tensor terms on top of the 
Skyrme interaction produce different effects. The experimental value of the Ml peak is 10.23 MeV [l^]. Since the 
theoretical values are 9.28 MeV (SLy5 without tensor), 12.31 MeV (SLy5 with tensor), 10.47 MeV (T44), it seems 
that the result obtained with T44 is preferable. 

We turn to the analysis of the results for ^°^Pb which are reported in Figs. [l]and[5l in the case of SLy5 plus the 
tensor force parameters of ^ and in the case of T44, respectively. In the unperturbed spectrum, the main role is 
played by the proton hii/2 — * hg/2 and neutron ii3/2 — >■ iii/2 configurations. They are clearly visible in the upper 
panels of the two figures. In the case of the SLy5 interaction, without tensor, the proton configuration lies at 5.85 
MeV while the neutron configuration is at 7.49 MeV. As already mentioned several times, in the case of SLy5 plus the 
the tensor parameters of [0] the spin-orbit splittings are increased when the tensor terms are taken into account: the 
energies become 6.45 MeV and 9.17 MeV, respectively, for the proton and neutron configurations we have mentioned. 

Within RPA, without tensor, we have found two peaks at 7.39 MeV and 9.14 MeV. The lowest peak is mainly 
composed by the proton hii/2 — > h9/2 configuration with an admixture of the neutron ii3/2 — > iii/2 configuration 
having different sign in its amplitude; the highest peak is mainly based on the neutron ii3/2 — > iii/2 configuration, 
with some admixture of the proton one, \iii/2 ^ ^'d/2i having the same sign in its amplitude. In other words, the 
lowest (highest) state has more IV (IS) character. This isospin character is not strongly pronounced because, with the 
values of T and U that have been employed, the non-diagonal matrix element which mixes the proton and neutron 
configurations is small. Since experimentally one finds that the lowest (highest) state has more IS (IV) character, 
there is a doublet inversion which is related to the values of the Landau parameters (discussed above) . We also notice 
that inversions of IS and IV spin doublet have also emerged from the self-consistent Gogny calculations of Ref. [l5| . 

With the inclusion of the tensor terms, the lowest and highest peak move, respectively, to 7.79 MeV and 10.57 
MeV. If we neglect the small mixing between the proton and neutron configurations, we can extract the values of 
the matrix elements of the residual tensor force from Eq. ^ separately for the proton and neutron states. We find, 
respectively, -0.20 MeV and -0.25 MeV. These values are smaller but have the same sign as in ^^Ca. 

Similarly to ^^Ca, the results in ^°*Pb obtained by using T44 seem better. In particular, we can notice that in this 
case one has the correct ordering in the spin doublet, the IS (IV) peak being the lowest (highest). In fact, with the 
values of T and U associated with this parameter set, the IS and IV characters of the states are more pronounced 
due to a larger value of the matrix element mixing neutrons and protons. The values of the energies are respectively 
6.12 MeV and 8.27 MeV and they compare better with the experimental findings (see e.g. [1^ and references therein) 
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which are 5.85 MeV and 7.30 MeV, than SLy5 without tensor (which gives 7.39 MeV and 9.14 MeV) or SLy5 with 
tensor (which gives 7.79 MeV and 10.57 MeV). 

IV. CONCLUSION AND PERSPECTIVES 

The tensor component of the bare nuclear force has been always of central interest for nuclear physics, but as far as 
effective interactions for many-nucleon systems are concerned, only recently we have achieved some first understanding 
of its role within self-consistent frameworks. In particular, in recent years many works have appeared in the literature 
which deal with the effect of tensor terms when they are added on top to, or fitted together with, Skyrme-type forces. 
These works have dealt with the impact of tensor terms on bulk properties like masses, or on single-particle states. 

Our work deals, for the first time, with the effects produced by these tensor terms on the multipole response of 
finite nuclei. We have attempted to disentangle the effects due to the modifications of the static mean field, and those 
due to the residual interaction, by analyzing the results of self-consistent HF plus RPA calculations performed in the 
case of quadrupole, octupole and magnetic dipole channels. 

The modifications of the static mean field had been already understood by means of our previous works. The 
dominant character (attractive or repulsive) of the residual interaction matrix elements has been extracted from the 
numerical calculations but also understood on the basis of a separable approximation for the tensor p-h force. Then, 
it has become evident that since these two effects are governed by different combinations of the parameters of the 
tensor force, the effects which are visible on the final results are ruled by a delicate interplay and are non-trivial. This 
is one of the main findings of the present work. 

The magnetic dipole states are more affected by the inclusion of the tensor terms at the mean field level, since 
unperturbed configurations correspond exactly to energy jumps between spin-orbit partners. The inclusion of the 
tensor residual interaction is demanded by self-consistency: for its matrix elements we have extracted typically values 
around few or several hundreds of keV, but these values depend of course on the values of the tensor force parameters 
T and U. 

Our work is exploratory and we have shown that the tensor force plays a role in RPA since the matrix elements are 
in general not negligible. We have considered two cases, namely the parameter set SLy5 plus the tensor parameters 
of [0] and the set T44. The latter seems preferable if one looks at the comparison with experiment in the magnetic 
dipole case, but not for some of the low-lying states in Ca isotopes. 

Actually, before learning more from the comparison with experiment, and fitting an ultimate Skyrme set including 
the tensor terms, one should ask the question whether the Skyrme ansatz is general enough for the spin-isospin channel. 
Some of the most modern, and most widely used, sets like the Lyon parameterizations have somewhat unsatisfactory 
values of the spin and spin-isospin Landau parameters and this anomaly should be cured. Improved Skyrme sets and 
Skyrme sets fitted with the tensor terms can be studied in detail, as far as their performance for excited states is 
concerned, by using the present RPA formalism. This is the main perspective opened by our present work. 

Acknowledgments 

C.L. acknowledges the support of the UniAMO fellowship provided by Fondazione Cariplo and Universita degli 
Studi which has allowed his stay in Milano, and the support of the National Science Foundation of China under Grant 
Nos. 10875150. C.L. and G.C. also acknowledge partial support from the Kavli Institute for Theoretical Physics 
China (KITPC). This work is partially supported by the Japanese Ministry of Education, Culture, Sports, Science 
and Technology by Grant-in- Aid for Scientific Research under the program number (C(2)) 20540277. 

APPENDIX A: CALCULATION OF THE P-H MATRIX ELEMENTS OF THE TENSOR INTERACTION 

We use the triplet-even and triplet-odd zero-range tensor terms which have been introduced originally by Skyrme 
[11] and read 

Ftensor = | |[(tTl • k' ) (^2 • k' ) - ^ (^i • a^)\l'']5{T^ - T^) + 5{v^ - T^)[{a^ ■ k) (^2 • k) - ^ (^i • (T2)k2]| 

+[/ ((ai • kV(ri - r2)(rT2 • k) - J(ai • a2)<5(ri - r2)[k' • k] j , (Al) 



where k = — acts on the left while k = ^'^„T^ acts on the right. Using the formula 
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(a . b)(c . d) = (f^^K^ + ^ " ^) + V (-l)-'>[a(^) « c(^)](^)[b(^) « d(^)]L^), (A2) 

q=-2 

we can re- write (jAip in such a way that the coupling of both the spatial and spin operators to rank 2 is evident, 



q= 2 J 

-C/| ^ (-ir-^[ai®a2](f)[k'0k]L')5(ri-r2)|. (A3) 



For RPA, we wish to calculate the coupled p-h matrix elements defined as 



ahcd 2^ 



^— {~iy^-'''^^''^"'H]ama,jc-m,\JM){3dma,3i>-mi\JM){amabmt,\^^^ (A4) 

We use the well-known Pandya relation to re-express the p-h matrix elements in terms of the particle-particle (pp) 



ones, 



yail = E j'\-iy^^'^^'' | } {Uajb)j'\v\u,u)j')^ (as) 



and change the jj coupling to the LS one, to obtain 



X 



{la L 4 1 ( Ic L Id 
Ja J 3b ){ 3c f Jd \{{LS)J'\V\{LS')J'). (A6) 
h S \ ] [ \ S' 



Then, 



{{LS)j'M'\V\{LS')j'M') = {-lf+''+'^' ^, ^ |(5||[ai®a2]('^||5')(L|| [6i®02]'"||i'>, (A7) 

where the operators Oi acts on the space part of the wavefunction. For the spin part, the reduced matrix element is 
simply 

[cTi ® (72]'^' \\S') = V3 X 3 X 5 X i X Ve X %/6 = 2\/5 (A8) 

9 

while the space part is quite lengthy to evaluate. After applications of the gradient formula and proper recouplings, 
one arrives at the final result, that is, 

{{LS)J'\V\{LS')J') ^ (-1)"+"'+"' { f f } X 

Q,A,A' 

{la h L ^ ^10 

Ic Id L' \ A'^A^ X LL' X 5 X - ^ (A9) 

A' a 2 J ^n=i 
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where 

X I I I ^ I (^crcll^a + ^ + 1 drr^Mcir)Mr)DlDt+^Mr), (AlO) 

w*^^^ has the same form as v'^^'' with the changes a < — > b and c < — > d, v^^^ has the same form as v'^^^ with the change 
a < — > c and an additional phase (—1)^ , and v'^^^ has the same form as w*-^-* with the changes a < — > d and c < — *■ b 
plus the additional phase (—1)^ Then 

i , IV A , ^ , 1 I 2 } f Ic + k la \ 

a la + i 



k,i=±l 

X 1 A' 

Ic k 

and v^^-^ has the same form as v^^^ with the changes a > 6 and c <— ^ cZ. Finally, 



Yla + i\\Yjlc + k){h\\YM I drr^Mr)Mr)DlMr)DtMr), (All) 



fe,i=±l 



2^ 2 J V 2 " 2; \ A 1 2 / \ + » 1 a 

(;a||>^a||/c + i)(/b||>^a||;d + fc> / drr^ 0, (r)06 (r) A\ (^c M A1 W , (A12) 



A 

Id + k 1 a 

and w'*^ has the same form as w^^^ with the changes c < — > a, d < — > 6 plus the additional phase (—1)^^'^ , w^^' has 
the same form as w'-^-' with the changes T — > U, c < — > a plus the additional phase (—1)^+'^, and u^^"-* has the same 
form as i;^^-' with the changes T — > U, b < — > d plus the additional phase (— 1)^"'*' . 
In all the formulas (jAlOp . (jAlip and (|A12p the following differential operator is used: 

dr r 



APPENDIX B: A SEPARABLE APPROXIMATION FOR THE TENSOR INTERACTION 

The tensor-even term in Eq. (|A3p can be written as 

Krnsor - |{[V5[f^i ® ^2]^') ® [k ® k']^^)]^ + V5[[ai ® (Ta]^') ® [k ® k](2)](")}5(ri - ra). (Bl) 

We assume that the contributions including gradients acting twice on the same wavefunction (~ Vi • Vi or ~ V2 • V2) 
are less important than the terms ~ Vi • V2 or ^ V2 • Vi, and we neglect them. With this approximation, and the 
multipole expansion of the (5-function, we can write one of the leading terms as 

or(tcrm 1) = I E ^^V)' ^ ^2]('' ^ [Vi V2](^)}(°)[rHl) ® Ye{2)fl (B2) 

e 

We apply standard angular momentum techniques to express this term of the interaction as a sum of separable terms, 
that is, terms made up with a product of two operators acting respectively only on particle 1 and particle 2. The 
result is 

,(term 1) = E (-)'+''+'AA'fc { \ ]] { y H] [Oa',.(1) ® 0.,,(2)](°), (B3) 
where the operator is 

dx,k{i) = W^ ® (V. ® Ye{i)Y^^]^''l (B4) 
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If we calculate the p-h matrix elements coupled to J [cf. the previous appendix, Eq. (|A4p ]. then k = J and 

^,''llp(term 1) = E (-)'+''+'AA'^ { i { y H] {p'\dyMh'){h\d,Mp)- (B5) 



We then analyze the quantity 

{p'\dx',k\h'){h\dx,k\p). (B6) 
Due to the time- reversal properties of the operator (|B4p , this is equal to 

- {p'\dyMh'){p\dx.k\h) (B7) 



(cf., e.g., p. 312-313 of 25]). To make an estimate of the sign of the diagonal matrix elements, we obtain from (|B5|) 
and 

^(J) ^fo._n-_^ ^ NA+A'+fcH'F f A A' 2 1 f 1 12 



e,x.x' 



^;i(terml) = -— E(-)'+''+'AA'fc } \ \ y x k (P|OA',fc|/»)(p|OA,fc|/i). (B8) 



We apply this latter equation for two typical cases (of interest for the low-lying spectroscopy of nuclei) , namely the 
natural parity case J^=2~^ and the unnatural parity case We obtain 

sign (even, term 1)) = -sign(r), (B9) 

sign (l^p'^'^^ (even, term 1)) = +sign(r). (BIO) 

Exactly the same kind of study can be performed for the tensor-odd part of the interaction, namely for the term 

^t°e'sor = C/{ V5[[ai ® a^r^ ® [k ® k](2)](0)}<5(ri - r^). (Bll) 

In this case, we assume that the leading term is ~ V'j^ • V2 (where the prime means that the gradient is acting at left). 
The equation corresponding to (jB2[) is 



y,°iL(term 1) = - j ^ ® ^2]^'^ ® [V'l ^-^f^Y'^Ml) r,(2)]W. (B12) 

By developing the same procedure outlined for the even term, we come to the conclusion that 

sign {V^hl (odd, term 1)) = +sign {U) , (B13) 

sign (yi,l2 (odd, term 1)) - +sign ([/) . (B14) 

We can use Eqs. (jB9|) . (|B10p . (|B13|) and (jB14p as a useful guideline. In fact, these equations correspond to direct 
(i.e., non antisymmetrized) matrix elements. The antisymmetrized interaction can be written as 

^tensor (1 " PAlPaPr) : (B15) 

where the position, spin and isospin exchange operators Pm, Pa and Pr have been introduced. Since the tensor 
interaction acts only among spin-triplet states, Pa=l- By definition, P]^j—+1 (-1) in the case of the even (odd) tensor 
term. Therefore, if we call for simplicity 

y(even) = ^^(even), (B16) 
y(odd) = l/^£(odd), (B17) 



it is simple to arrive at 

yihPp = Qv^(even) + ^y(odd) ) IS part 

-iv^(even) + iv^(odd) ) (n • T2) IV part. (B18) 
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The sign of the residual interaction can be understood by combining this latter equation with the previous ones. 
This is at least true when the different terms add up coherently; in case of cancellations, of course, the qualitative 
arguments of this Appendix may not be valid. 
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TABLE I: Parameters of the effective interactions used in the HF plus RPA calculations. The set SLy5 has been introduced 
in [3], and it is here supplemented by the tensor force parameters T and U which have been introduced in The set T44 is 
taken from [8|. The notation is the same as in [H for the Skyrme parameters, while T and U are defined as in f3|. 





SLy5+(T,U) 


T44 


to [MeV fm^] 


-2484.88 


-2485.670 


ti [MeV fm'^] 


483.13 


494.477 


t2 [MeV fm'^] 


-549.40 


-337.961 


^3 [MeV fm^+S"^] 


13763.0 


13794.7 


Xo 


0.778 


0.721557 


Xl 


-0.328 


-0.661848 


X2 


-1.000 


-0.803184 


X3 


1.267 


1.175908 


a 


1/6 


1/6 


VKo [MeV fm'^] 


126.0 


161.367 


T [MeV fm'"^] 


888.0 


520.983 


U [MeV fm^] 


-408.0 


21.522 



TABLE IL Properties of the low-lying quadrupole states in ''^Ca and ^°*Pb. Energies and B(E2) values from the RPA 
calculations are compared with experimental data from Ref. 





Energy [MeV] 


B(E2) [e^fm"] 


«Ca 


SLy5 (no tensor) 


3.05 


56.5 


SLy5 (with tensor) 


3.88 


52.9 


T44 (no tensor) 


3.35 


53.2 


T44 (with tensor) 


3.17 


51.8 


Exp. 


3.83 


95.0 


208 pj^ 


SLy5 (no tensor) 


4.89 


2.94 10^ 


SLy5 (with tensor) 


5.06 


3.07 10^ 


T44 (no tensor) 


5.04 


2.86 10^ 


T44 (with tensor) 


4.49 


2.70 10^ 


Exp. 


4.07 


2.97 10^ 
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TABLE III: The same as Table HIl for the case of the low-lying octupole states. In this case the nuclei are «Ca, 4«Ca and ^ospb. 
The experimental data are from Ref. p^ . 

Energy [MeV] B(E3) [e^fm^] 

*°Ca 

SLy5 (no tensor) 3.78 1.32 10* 

SLy5 (with tensor) 3.02 0.91 10* 

T44 (no tensor) 1.36 1.23 10* 

T44 (with tensor) 1.22 1.03 10* 

Exp^ 3/74 2.04 10* 

*»C^ 

SLy5 (no tensor) 4.78 0.54 10* 
SLy5 (with tensor) 6.16 0.66 10* 
T44 (no tensor) 3.44 0.42 10* 
T44 (with tensor) 4.93 0.34 10* 
Exp^ 4^51 1.00 10* 

SLy5 (no tensor) 3.51 6.90 10^ 

SLy5 (with tensor) 3.49 6.16 10^ 

T44 (no tensor) 3.49 5.87 lO'' 

T44 (with tensor) 3.19 6.29 lO'^ 

Exp. 2.61 (5.30±0.30)10'^ 
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FIG. 1: Unperturbed (left panels) and RPA (right panels) strength functions associated with the isoscalar quadrupole operator, 
in the case of ^"Ca, **Ca and ^"^Pb. We display results both with and without the inclusion of the tensor force in the case of 
the force SLy5 plus the T and U parameters of 6] (first column of Table H]). The discrete RPA peaks have been smeared out 
by using Lorentzian functions which have 1 MeV width. 
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FIG. 2; The same as Fig. [T]in the case of the octupole response. 
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FIG. 3: Unperturbed and RPA strength functions in ^®Ca associated with the isoscalar and isovector Ml operators defined by 
Eqs. Q and We display results both with and without the inclusion of the tensor force, in the case of SLy5 plus the T 
and U parameters of Q (left panels) and in the case of the T44 parameter set (right panels). The discrete RPA peaks have 
been smeared out by using Lorentzian functions which have 1 MeV width. 
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FIG. 4: Unperturbed (upper panel) and RPA (middle and lower panels) strength functions in ^"*Pb associated with the isoscalar 
and isovector Ml operators defined by Eqs. ([Sj and Q. We display results both with and without the inclusion of the tensor 
terms, in the case of the force SLy5 plus the T and U parameters of @| (first column of Table The discrete RPA peaks have 
been smeared out by using Lorentzian functions which have 1 MeV width. 
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